Quasars as bubbles of dark matter: evidence for axion and tachyon matter in the 

Universe 



Anatoly A. Svidzinsky 
Department of Physics, Institute for Quantum Studies, Texas A&M University, TX 77843-4242 

(Dated: February 2, 2008) 

Growing amount of data show evidence for statistical and apparent physical association between 
low-redshift galaxies and high-redshift quasi-stellar objects, suggesting noncosmological origin of 
their redshift and failure of classical quasar explanation. Here we find an analytical solution of 
Einstein equations describing bubbles made from axions with periodic interaction potential. Such 
particles are currently considered as one of the leading dark matter candidate. The bubble interior 
has equal gravitational potential and, hence, photons emitted from the interior possess identical 
gravitational redshift. The redshift depends on the bubble mass and can have any value between 
zero and infinity. Quantum pressure supports the bubble against collapse and yields states stable 
on the scale of the Universe age. Our results explain the observed quantization of quasar redshift 
and suggest that intrinsically faint point-like quasars associated with nearby galaxies (a few % of 
known objects) are axionic bubbles with masses 1O S -1O 9 M and radii 10 3 -10 4 i? Q . They are born 
in active galaxies and ejected into surrounding space. Properties of such quasars unambiguously 
indicate presence of axion dark matter in the Universe and yield the axion mass m — 0.4 — 3 meV, 
which fits in the open axion mass window constrained by astrophysical and cosmological arguments. 

We also found that tachyons, another dark matter candidate, can form objects with galactic size, 
negligible mass and any value of the gravitational redshift. Such finding combined with quasar 
observations suggests that bright quasars 3C 48, 3C 273 and 3C 279 are nuclei of forming nearby 
small galaxies embedded into tachyonic clots and possess pure gravitational redshift. If the bright 
quasars later evolve into small companion galaxies, then their dark galactic halos, observed by 
rotation curves, are probably remnants of the tachyon matter. 



I. INTRODUCTION 

Since their discovery in 1960's quasars (QSOs) be- 
came one of the most mysterious objects in the Universe. 
All observed QSOs share several remarkable traits. First 
they are strong sources of radio waves and their radia- 
tion spectrum is not a Planck spectrum of thermal radi- 
ation as in ordinary stars. Secondly, all QSOs produce 
radiation with significant redshift, some of them possess 
redshift z > 6. Third, the radiation is characterized by 
strong emission lines, while ordinary stars produce light 
with strong absorption lines. Finally, all QSOs have small 
volumes by galactic standards. One of the most puzzling 
features is that if the redshift occurs due to large speed of 
QSOs relative to Earth, they must be cosmological dis- 
tances away and extremely luminous - their luminosity 
should lie between 10 and 100 of the total optical lumi- 
nosity of the brightest galaxies. 

The present-day understanding of QSOs is based on 
the paradigm that redshift of all QSOs has cosmological 
origin (the so called cosmological hypothesis). Implied 
huge quasar's luminosity is the main consequence of the 
cosmological hypothesis. Based on this assumption one 
can find some similarity between QSOs and active galac- 
tic nuclei: both are intense energy producers, both pos- 
sess strong emission from radio to X-ray bands, some- 
times both have jets issuing from the central region. As 
a result, the present-day conventional view regards QSOs 
and active galactic nuclei as having the same nature. 

The classical (generally accepted) picture of QSOs as 
active galactic nuclei (AGN) is based on the hypothesis 



proposed by Sandage that all QSOs lie in the nuclei 
of galaxies. The weaker of them are the classical Seyfert 
1 nucleus, where the galaxy of stars surrounding the nu- 
cleus is clearly seen. In the unified scheme the QSOs are 
high luminosity Seyfert 1 nuclei which contain a compact 
nuclear source ionizing a broad line region, surrounded 
by an optically thick torus of dust. Depending on the ori- 
entation of this torus with respect to the line of sight, the 
central object is seen or hidden. When it is hidden, we 
see only the narrow, extended emission line region; the 
galaxy is a Seyfert 2. When it is seen, we see both broad 
and narrow emission lines. Most radio loud sources have 
a double lobe structure which is powered by a relativistic 
jet. When the angle between the jet axis and the line of 
sight is small the object appears as a blazar. 

However, despite of great efforts, the nature of QSOs 
remains unclear so far. The point is that the unifying 
scheme in its present form contradicts to some obser- 
vations which can not be ignored. Problems with the 
current quasar's paradigm are discussed in many publi- 
cations and here we mention the main of them (for details 
see, e.g.. 

H USE 00)- One of the consequences of 
the unified scheme is that host galaxies of all QSOs must 
be present. However, observations with Hubble Space 
Telescope [1 E3 have shown that for some QSOs there 
are no host galaxies, while for others there is evidence 
for faint galaxies or irregular systems with the same red- 
shift. Another problem is distances to quasars. Some 
QSOs and the compact galaxies they evolve into are re- 
solved and many low redshift QSOs lie close to g alaxies 
with approximately the same redshifts PEIII. This, 
as usually believed, suggests on cosmological distances 
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to the objects. On the other hand, there is a strong 
evidence that many low redshift (nearby) galaxies and 
high redshift QSOs are physically associated and, hence, 
these QSOs are no further away than the close galax- 
ies and must have redshifts noncosm olog ical in origin 
(see, e.g., @ H H Q E E [H E M, HI and ref- 
erences therein) . Observations suggest that such quasars 
are ejected from active galaxies or in the process of ejec- 
tion from the galactic nucleus. 

A systematic search for quasar-galaxy association was 
made by Burbidge et al. They found over 500 

close pairs of QSOs and galaxies of which only 28 QSOs 
were associated with 42 galaxies with approximately 
equal redshifts (this is what the cosmological hypoth- 
esis predicts), while for the reminder 92% objects the 
quasar redshifts considerably exceeded the galaxy red- 
shifts (zg a i = 0.0001 — 0.01). Another evidence is the 
observation of three QSOs with different large redshifts 
(z = 0.60, 1.40, 1.94) near a spiral arm of the galaxy 
NGC 1073 by Arp and Sulentic 0. The galaxy NGC 
1073 has much smaller redshift z = 0.004 and approxi- 
mately 16 Mpc away (for H — 75km/s-Mpc). There are 
dozens of similar evidences, some of them are discussed in 

HBEEmEiim. Itis im p° rtant to note that there 

is a significant excess of galaxies only around flat spec- 
trum radio- loud QSOs (blazars), which constitute about 
10% of the known quasar population. In contrary, there 
is a marginal defect of galaxies around optically selected 

QSOs nEm. 

One of the most fundamental quasar property is the 
redshift quantization. Based on observations, Karlsson 
|24j has noted division of quasars into two groups with 
different redshift properties and concluded the following. 
If we select QSOs associated with most nearby (distance 
d < 50 — 100 Mpc), galaxies then their redshift is close 
to certain values (quantized), as shown in Fig. 0]in next 
section. Meanwhile, in QSO samples associated with dis- 
tant, d > 50 — 100 Mpc, galaxies no periodicity in intrin- 
sic redshift is observed. Such a division is supported by 
later studies of QSOs associated with most near by g alax- 
ies where the quantization was confirmed [2^ and 
distant (0.01 < z ga i < 0.3) galaxies for which absence 
of any periodicity was claimed |27| . The property sug- 
gests existence of intrinsically faint (optical luminosity 
L = 10 5 — 10 7 i©) QSO subgroup with quantized non- 
cosmological redshift. Being intrinsically faint, such ob- 
jects are not detected from large distances which yields 
disappearance of redshift quantization in distant QSO 
samples. 

The observational evidences indicate that the total 
quasar population must be divided into several groups 
having different nature but similar radiation mechanisms. 
Quasars of the first group are intrinsically faint point- 
like objects in the optical band and have structure un- 
resolved in telescopes. Such objects are born in nearby 
active galaxies and ejected into surrounding space. Their 
redshift is quantized and noncosmological (mostly grav- 
itational). Here we show that bubbles of dark matter, 
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FIG. 1: Structure of intrinsically faint quasars 

pictured in Fig. ^ with masses about 10 8 — 10 9 M© and 
radii 10 3 — lO 4 i?0 can explain the intrinsically faint point- 
like quasars. The bubble is supported against collapse by 
quantum pressure and, in principle, can have decay time 
larger then the age of the Universe. Hypothetical ax- 
ions, one of the leading dark matter candidate, fit well 
into this picture and can account for the redshift quanti- 
zation. Usual baryonic matter falls into the bubble inte- 
rior, heated by the release of the gravitational energy and 
produce electromagnetic radiation. The bubble can also 
harbor a black hole at the center which makes its struc- 
ture similar to AGN. The amount of baryonic matter 
trapped in the bubble we assume to be small compared 
to the bubble mass. Bubble radiation spectrum must 
be similar to AGN because the radiation mechanism is 
the same. Photons emitted anywhere inside the bubble 
interior posses identical gravitational redshift and freely 
propagate into surrounding space because dark matter is 
transparent for electromagnetic waves. 

Quasars of the second group are intrinsically bright ob- 
jects with typical optical luminosity L = 10 8 — 10 w Lq 
(a fraction of the objects can be more luminous). Some 
of them are resolved in telescopes. Their redshift pos- 
sesses both the intrinsic and cosmological components, 
but the intrinsic redshift is not quantized. In Section 4 
we show that not only small on galactic scale objects can 
posses large gravitational redshift. Droplets of tachyon 
dark matter with the mass much smaller then galactic 
mass can produce any value of the gravitational redshift 
on kpc scales. Baryonic matter embedded in such re- 
gions behaves as being possessed of an intrinsic redshift. 
Tachyonic clots can explain the nature of intrinsically 
bright QSOs which possibly have noncosmological red- 
shift component but occupy large space volumes 1 — 100 
pc in size. The classical model of quasars, as active galac- 
tic nuclei with pure cosmological redshift, is a limiting 
case of tachyonic quasars when the gravitational redshift 
component is negligible compared to the cosmological 
contribution. Intrinsically bright quasars are observed 
from substantially greater distances then the faint ax- 
ionic bubbles and, as a consequence, dominate in the 
known quasar population. We discuss this issue in the 
last section. 

Luminous active galactic nuclei with pure cosmological 
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redshift belong to the third group of quasars. However, 
as we mentioned, there is no sharp distinction between 
such objects and the tachyonic quasars. 

Being intrinsically faint, the axionic bubbles constitute 
only a few % of the observed QSOs population. Intrin- 
sically bright QSOs are possibly a heterogeneous family. 
Probably blazars, about 10% of the known QSOs, have 
the tachyonic nature and their redshift possesses substan- 
tial gravitational component. First quasars discovered 
belong to this category. At the same time, maybe most of 
the radio-quiet (steep radio spectra) QSOs (about 90% of 
the known QSOs population) are luminous active galac- 
tic nuclei at the cosmological distances indicated by their 
redshifts. However, the last question requires further 
analysis. 

For the first time our theory was presented at the Fifth 
International Heidelberg Conference on Dark Matter in 
Astro and Particle Physics (3-9 October, 2004) [2^1; in 
this paper we provide all the details. 



II. DARK MATTER BUBBLES 

There is abundant evidence that the mass of the Uni- 
verse is dominated by dark matter of unknown form 
|29j . Particle dark matter, i.e., one or more relic par- 
ticle species from the big bang, is strongly suggested 
as the dominant component of matter in the Universe 
|30|. For many years weakly interacting massive parti- 
cles (WIMPs), such as, for instance, the lightest neu- 
tralino in the minimal supcrsymmetric standard model, 
were the leading dark matter candidate. However, there 
is still no evidence in favor of WIMP, either from bolome- 
ter experiments designed for direct detection or from the 
observation of cosmic rays. Moreover, many recent sim- 
ulations of structure formation in the Universe suggested 
that any dark matter component modelized as a gas of 
free particles, such as WIMPs, result in cuspy density 
profiles at galactic centers, while observation of rotation 
curves indicate a smooth core density [3l| . Such a contro- 
versy suggests that the dark matter halo probably con- 
sists of a Bose condensate of light particles which be- 
haves as a classical scalar field, coherent on the scale of 
10 kpc. Self gravitating Bose condensate, governed by 
the Klein-Gordon and Einstein equations, could account 
for the dark matter distribution inside galaxies (see, e.g., 
[UlillllEIIfil and references therein). The field acts 
as an effective cosmological constant (dark energy) be- 
fore relaxing into a condensate of nonrelativistic bosons 
[37l l38[ . In order to obtain the Bose condensate halo 
on the scale of 10 kpc one should consider ultralight 
nonintcracting particles with masses m ~ 10 -23 eV or 
heavier particles with self-coupling |33, |39(- From the 
rotation curve of the Andromeda Galaxy an estimate of 
the mean mass density of the dark matter composing the 
30 kpc luminous core is p ~ 2 x 10~ 24 kg/m 3 40]. Par- 
ticle concentration p/m determines the temperature of 
Bose condensation transition T c . The condition that T c 



is greater than the temperature of the cosmic microwave 
background radiation imposes an upper limit on the par- 
ticle mass m < 1 eV. So, in order to explain the dark 
matter halo by a Bose condensate of scalar bosons the 
particle mass must be in the range m = 10~ 23 — 1 eV. 

Here we study massive real scalar field ip with periodic 
interaction potential 



V(<p) = V [l-aos(<p/f)], 



(1) 



where Vq > 0. This potential is quite general and de- 
rived in quantum filed theory in connection with pseudo 
Nambu-Goldstone bosons (PNGBs) j^l El EH In all 
such models, the key ingredients are the scales of global 
symmetry breaking / and explicit symmetry breaking 
(lo) 1 / 4 . From the viewpoint of quantum field theory, 
PNGBs are the only way to have naturally ultralow 
mass, spin particles. One of the example of a light 
hypothetical PNGB is the axion which possess extraor- 
dinarily feeble couplings to matter and radiation and is 
well-motivated dark matter candidate |3*flj . Axion arises 
from the Pcccci-Quinn solution to the strong CP problem 
|43j . If the axion exists, astrophysical and cosmological 
arguments constrain its mass m to be in the range of 
10~ 6 — 3 x 10~ 3 eV and the global symmetry-breaking 
scale to lie in a narrow window 



Tfl-K fir 

2m 



= 2 x 10 9 - 6 x 10 12 GeV, 



(2) 



where m w = 135 MeV is the neutral pion mass and 
/tt = 93 MeV its decay constant [30l |44| . Axions in this 
mass range could provide much or all of the cold dark 
matter in the Universe. Interaction of axions with QCD 
instantons generates the axion mass and cosine interac- 
tion potential [45J. Another example of ultralight PNGB 
is a hypothetical scalar field that arises when the global 
symmetry is broken on a Planck scale / ~ 10 18 GeV and 
explicit breaking scale is comparable to light neutrino 
mass (Vb) 1 / 4 ~ 10~ 3 eV. Such a field, which acquires a 
mass m ~ y/Vo/ f ~ 10~ 24 eV, would currently dominate 
the energy density of the Universe [37j . 

A self-gravitating real scalar field tp in general relativity 
is described by the action 46] (in this section we use 
natural units for which h = c = 1) 



1= d 4 



1 



16ttG 



:Ri- 



(3) 

where G is the gravitational constant, J is the deter- 
minant of the metric tensor g^ v and Ri is the Riemann 
curvature. We consider spherically symmetric system, so 
the metric can be written in the form 



ds 2 



-N 2 dt 2 + g 2 dr 2 



r 2 dtt 2 , 



(4) 



where g, the radial metric, and N, the lapse, are functions 
of t and r with r being the circumferential radius. We 
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introduce dimensionless coordinates and define the unit 
of distance, time and ip as 



with the following boundary conditions 



ro 



h 

mc 



in = 



1 



Po = 



(•5) 



where c is the speed of light, m = yVo/ f is the mass 
of the scalar particle. In dimensionless units the coupled 
Klein-Gordon and Einstein equations describing dynam- 
ics of the field cp and the metric are |4(j 



g(0) = g{oo) = N(oo) = 1, 



g'(0) = AT'(O) = </(0) = 0, V(p(oc)) = 0. 

Outside the bubble <p = and Eqs. {TUJ-fEJl lead to the 
Schwazschild solution for a spherically symmetric prob- 
lem: 



f 

N 2 



Nip 
~N3 



p' fg 2 + l 



dV 

dp 



rip p 
~N 2 ~ 



2rg 2 V 



0, 



(6) 
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1 - 2M/r ' 



N 2 = l-— , (13) 



where M is the dimensionless bubble's mass in units of 
mlJm. 

Eqs. I|1U|) - (|12[1 can be obtained as an extremum condi- 
tion of the energy functional 



2 



- 1 



(7) E[N,g,ip] 



dr — | r 2 ( — 
9 I V 2 



/2 



:(5~1) 2 + 



1 - <? 2 



(S) 



where an overdot denotes d/dt, a prime denotes d/dr 



V 



:[1 — cos(a</j)], a 



(9) 



is the dimensionless potential and the coupling param- 
eter respectively, m p i = y/hc/G = 1.2 x 10 19 GeV is 
the Planck mass. The interaction potential V has de- 
generate minima at ip — 2irn/a, where n is an integer 
number. It is tempting to search for time-independent 
solutions of Eqs. JBJ-JSJ). However, the pseudovirial the- 
orem of Rosen _47] implies that no such solution is pos- 
sible in the Newtonian limit, and in the strong-field case 
it has been shown numerically that no nonsingular solu- 
tion exists . Here we show that in the limit of strong 
nonlinearity (a S> 1) an approximate time- independent 
solution does exist and it describes a spherical bubble 
with surface width much smaller then its radius R. The 
bubble surface is an interface between two degenerate 
vacuum states with p = 2im/a (r < R) and <p = 
(r>R). 

We look for static solutions, then Eqs. ©-(JSJ) reduce 

to 



dV 



g 2 + l , \ p" 

y —— - 2rg 2 V + ^ - — = 0, 
r J g z ops 



(10) 



"(<? - 1) 



N~ 



(14) 



Variation of this functional with respect to <p gives the 
Klein-Gordon equation l|l(J|) , while variation with respect 
to the metric functions g and N yields the Einstein equa- 
tions l[TT|l. JT2J respectively. Using Eqs. (jT2J), C3} the 
total bubble energy reduces to 



E 



dr 



1 

1 - - 

9 



N 



= lim r 

r — >oo 



1 



l--j=M, 
(15) 

here we applied the Schwazschild solution l|13|) valid out- 
side the bubble. Eq. (|15fl demonstrates Einstein equiv- 
alence principle between the mass and the energy. Eq. 
(|12fl leads also to another expression for the bubble mass: 



M 



dr 



1 

V 2 



drr 



2V 



(16) 

which shows that the scalar field gradient and the poten- 
tial V are sources of the bubble energy. 

Let us assume that a ^> 1 and the radius of the bubble 
is R 3> a. Then inside the bubble, including its surface 
region, one can omit terms with 1/r in Eqs. i|l(J[l - l|12|) 
and take r ~ R (thickness of the bubble surface is much 
smaller than its radius), we obtain 



2R(p'V 



dV 
dp 



0, 



(17) 



N>=» 
2 



9 - 1 



+ r (ip /2 - 2g 2 V) 



(11) 



(18) 



■9 = o 



1-9 2 



2g 2 V) 



(12) 



g> = 9 4 (p' 2 + 2g 2 V) 



(19) 
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Eqs. (|17|l - i|19|) can be solved analytically. Their first 
integral is 



Further, Eqs. 



-/V = const, 

tp' 2 = 2<?V, 

g' = Rgtp' 2 . 
and (0 yield 

■<p(o) 

VzVdtp. 



~ = 1-R 



(20) 
(21) 
(22) 

(23) 



We assume that <p(0) = 2irn/a, where n — 1,2,3,... 
is the number of kinks at the bubble surface, and ip(r) 
monotonically decreases with r. Outside the bubble tp = 
0. Substitute (J23J into J23J leads to 



V = 



1-Rf 



v(o) 

V 



(24) 



Further integration can be made for a particular choice 
of the scalar field potential. For V(ip) given by Eq. © 
the final solution is 



AR 

— In | sin(ay>/2)|- 



^(2m-l) 



arctanh[cos(a(,c/2)] 



= sign[sin(a</?/2)](r - R m ), 
tp G [2n(n — m+ l)/a. 2n(n — m) / a] 



(25) 



where R m is a position of the mth kink, m = 1, 2, . . . , n. 
When the coordinate r passes through the point R m 
the scalar field tp(r) changes from 2ir(n — m + l)/a to 
2ir(n — m)/a (see Fig. 2a). Eq. lj2l5|) yields the following 
expression for g as a function of tp inside the bubble: 



i = l-l#[2^-l 
g a z 



■ cos(q!(/j/2)]. 



(26) 



Outside the bubble tp = 0, m = n and X/g = 1 — 
8nR/a 2 . The solution is valid if 1/g > 0, that is 
R < -R m ax = a 2 /8n. Match of the inner solution 12tj|) 
with the Schwazschild solution (|13|) outside the bubble 
determines the bubble mass: 



M 



8nR 2 32n 2 R 3 



(27) 



The mass depends on the bubble radius R and changes 
between zero for R = and a maximum value M max = 
a 2 /16n for R — R ms , x . 

The mass-radius relation l|27|) was derived for the co- 
sine interaction potential ©. However, using Eq. 1)23(1 . 
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one can obtain similar relation for any periodic potential 
Vbp) 



M = AnnuR 2 - 8?r 2 n 2 w 2 i? 3 , 



(28) 



where u is the surface energy density (the energy per unit 
area) determined by an integral over one potential period 



— / V2Vdtp. 

47T ' 



(29) 



For the cosine potential JSJl u = 2/ ira 2 . 

Redshift of the bubble interior z — l/N — 1 can be 
found by matching the inner solution TV = const with 
the Schwazschild solution (|13|) outside the bubble. As a 
result, we obtain that everywhere inside the bubble the 
redshift is the same and given by 



1 



yjl - 2M/R 



1 



1 



1 — AnnuR 



1. 



(30) 



The internal redshift monotonically increases from zero 
to infinity when the bubble radius R changes from zero 
to -R max - Fig. 2b shows redshift of space as a function of 
the distance r to the bubble center. The redshift is con- 
stant in the bubble interior and monotonically decreases 
outside the bubble. 

In Fig. |31 we plot the bubble gravitational redshift as 
a function of its mass for different kink number n. The 
redshift monotonically increases with the mass and goes 
to infinity when the mass approaches the maximum value 



a /IQn. 



A. Redshift quantization 

Let us make rescaling M — > M/Attu, R — > R/Airu, then 
Eqs. lO, © become 



M = nR 2 - n 2 R 3 /2, 
1 1. 



1 - nR 



(31) 
(32) 
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FIG. 2: (a) Scalar field p> as a function of distance r to 
the bubble center for bubbles with equal radius and differ- 
ent quantum numbers n = 1, 2, 3. The unit of length is 
h/mc. Note, we plot the field ip only in the vicinity of the 
bubble surface where it undergoes variation, (b) Redshift z 
of space as a function of distance r to the bubble center for 
bubbles shown in Fig. 2a. The redshift is constant in the 
bubble interior. 




FIG. 3: Redshift of a bubble as a function of mass for different 
"quantum" numbers n =1, 2, 3, 4, 5, 6. 



For a given M the redshift depends on the number of 
kinks n, which implies that the redshift is quantized. 

In early QSO samples involving about 600 radio- 
emitting QSOs associated with bright nearby spiral 
galaxies (z ga i < 0.007), Karlsson showed that the redshift 
distribution has a periodicity log(l+z K +i) — log(l+z n ) = 
0.089, where n = 0,1,2, ... and z = 0.061 JMj3. It 
has been later confirmed by other groups |25l l50| . In 
a recent paper, Burbidge and Napier |2(| tested for the 
occurrence of this periodicity in new QSO samples and 
found it to be present at a high confidence level. The 
peaks were found at z « 0.30, 0.60, 0.96, 1.41 and 1.96 
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R, in a. 2 h/mc 


z 


1 


0.0329 


0.357 


2 


0.0241 


0.629 


3 


0.0204 


0.96 


4 


0.0182 


1.40 


5 


0.0168 


2.06 


6 


0.0159 


3.24 


7 


0.0153 


6.11 


8 


0.0151 


26.6 



TABLE I: Redshift of the bubble interior z and its radius R 
for M = 0.00752a 2 m^/m and different kink numbers n. 

in agreement with Karlsson's empirical formula. The for- 
mula also includes the peak at zq = 0.061, however, this 
peak does not occur for quasars, but for morphologically 
related objects. 

The redshift periodicity is observed only in QSO sam- 
ples satisfying certain selection criteria, in particular, the 
galaxies which are assumed to be paired to the QSOs 
must be most nearby spirals p4L lRlj | . This implies that 
redshift quantization is a property of intrinsically faint 
QSOs which are not detected from large, > 50 — 100 
Mpc, distances. The observed peaks are narrow, imply- 
ing that the spread in cosmological redshift and Dopplcr 
redshift due to random motion must both be small. If 
they were not, the intrinsic peaks would be washed out, 
as is easily seen from the combined redshift formula 

1 + 2= (1 + z g )(l + z c )(l + z d ), 

where z, z g , z c and Zd denote the total, gravitational, 
cosmological and Doppler redshift respectively. The mor- 
phology of companion galaxies is probably also important 

If we consider a sample of quasars born in the same 
type of galaxies it is naturally to expect that such objects 
would have approximately equal masses because their for- 
mation mechanism must be similar. Such phenomenon is 
well known for type la supernovae or neutron stars: prac- 
tically all measured neutron star masses cluster around 
the value of IAMq with only a few percent deviation 
[52|. If we assume that dark matter bubbles are born 
with equal masses then their redshift must be quantized. 
For a given M Eqs. (|31[) . (|32|l have a set of solutions for 
R and z corresponding to different values of the "quan- 
tum" number n. For M = 0.0601 (in dimension units 
M = 0.00752a 2 m 2 , 1 /m) Eqs. © have solutions for 

n = 1, 2, . . . , 8, they are given in Table 1. 

In Fig. 2] we plot the most recent histogram of the 
redshift distribution from Ref. [5l| in which five peaks 
are clearly seen. The solid lines show the redshifts from 
our Table 1, they match well the observed peaks. The 
agreement is remarkable because the theory has only one 
free parameter, the bubble mass M. Such coincidence 
strongly suggests that some fraction of quasar popula- 
tion is dark matter bubbles composed of scalar particles 
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FIG. 4: Histogram of the redshift distribution of QSOs close 
to bright active spiral galaxies or multiple QSOs with small 
angular separation from Ref. |Hl|. The solid lines represent 
position of the peaks from Table 1. 



with periodic interaction potential. One should mention 
an alternative possibility of quasar evolution. Bubbles 
can be originally born with the same mass and number 
of kinks n = 5 that corresponds to the 5th peak. Dur- 
ing evolution the kinks tunnel to the bubble center and 
quasars sequentially decay into states with smaller n but 
the same mass. 

If a bubble is made of axions with m = 10" 4 -3x 10" 3 
eV and / w f^m^/2m = 2 x 10 9 - 6 x 10 10 GeV, 
then, according to Eq. ©, a = 5.6 x 10 7 - 1.7 x 10°. 
Then a bubble with the internal redshift z — 0.36 and 
n = 1 would have the mass M — 0.00752a 2 m 2 1 /m = 
3 x 10 7 - 10 9 M o , surface width ~ h/mc = 0.07 - 2mm, 
surface mass density 4 x 10 12 — 10 14 kg/m 2 and the radius 
R = 0.0329a 2 n/mc = 3x 10 2 -10 4 i? Q . Such radius range 
agrees with the size of the emission region expected for 
the intrinsically faint point-like quasars associated with 
nearby galaxies. Indeed, analysis of the emission lines 
intensity of the bright quasar 3C 48 shows that if it was 
located a distance of 1000 Mpc then the size of the broad 
emission-line region should be about 0.1 pc 53] . The 
quasars we consider as candidates into dark matter bub- 
bles have the apparent visual magnitude 18-21 0, ITsj 
which yields the brightness 1-2 orders smaller that 16th 
magnitude quasar 3C 48. If they are located 1 — 10 Mpc 
away and radiating gas has the same parameters as for 
3C 48 their size must be of the order of 10 4 i?© (we as- 
sume that the brightness scales as R 3 /d 2 , where d is the 
distance to the object). 

An empirical relation between the size of the broad- line 
region R and luminosity L of Seyfert 1 galaxies, R oc L 1 ^ 2 
|54j yields a similar answer. Indeed, for Seyfert 1 galaxies 
R ~ 10 — 100 light days. The luminosity of the point-like 
quasars is 5-6 orders smaller, which for the quasars leads 
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FIG. 5: Redshift of the bubble interior z as a function of time 
t during collapse of initialy static one kink bubble with the 
radius Ro = 0.25a 2 /8 and a = 1000. The unit of time is 
h/mc 2 . 



to R ~ 10 4 i? ( 



III. STABILITY 

The obtained analytical solution of the stationary 
equations I|l() |l -p2 |l is approximate because the terms of 
the order of 1/a 2 were omitted. For a 3> 1 such terms are 
very small, however, they play a crucial role in the bub- 
ble stability. Our numerical calculations show that the 
complete system of stationary equations l)10 [l -l)12 fl has no 
solution which indicates a possible instability. To study 
the unstable mode we solve numerically the evolution 
equations for the scalar field ©-JSJ with the initial con- 
dition given by the approximate analytical solution 12511 , 
(|26[) : the result is shown in Figs. |S] and HJ1 Under the 
influence of surface tension and gravitational attraction 
the initially static bubble starts to collapse very slowly 
with the acceleration a = R ~ —c 2 /R (the corresponding 
unstable mode has an imaginary frequency u> ~ ic/R). 
The contraction is accompanied by radiation of scalar 
particles from the bubble surface which propagate into 
surrounding space. One can describe the instability as a 
tendency to reduce the surface area of a hollow bubble 
with the surface tension a = u. 

The dynamics of relativistic bubbles has been investi- 
gated in connection with the bag model of hadrons and 
phase transitions in the early Universe [55ll56ll57ll58ll59| . 
The study has mostly focused on the thin-wall approxi- 
mation where the bubble surface is specified by effective 
parameters such as the surface energy density u and the 
surface tension a. The approximation omits the particle 
emission and describes the bubble dynamics by a single 
differential equation for the radius R. In particular, when 
R = the acceleration is given by [5(j 
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FIG. 6: Radius of the bubble R as a function of time t dur- 
ing collapse of initialy static bubble with the radius Rq — 
0.25q 2 /8 and a = 1000. The unit of time is Ti/mc 2 . 
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(33) 



where N = y/l - 2M/R. The first term on the right- 
hand side of Eq. (|33|l is the Laplace pressure, while the 
second term is the pressure caused by gravitational at- 
traction. At zero temperature u — a. Our numerical 
result for the surface acceleration (valid at T = 0) agrees 
well quantitatively with Eq. Ij33(l obtained in the thin- 
wall approximation. This indicates that energy dissipa- 
tion due to emission of scalar particles is negligibly small 
and does not affect the bubble contraction. 

When the bubble size approaches the gravitational ra- 
dius R g = 2M a distant observer sees that shrinking bub- 
ble starts to slow down and the shrinking speed decreases 
to zero. Fig. [7| shows solution of the evolution equations 
©-© in the regime when such a behavior is achieved. 
The bubble radius never crosses the event horizon; in the 
reference frame of the distant observer the radius R(t) 
only asymptotically approaches R g . At large t the bub- 
ble radius R(t) is well approximated by an exponentially 
decaying function shown by dashed line in Fig. [7] The 
effect is caused by time dilation. Motion of a particle 
falling into a black hole with a gravitational radius R g 
is a classical example of such effect. A distant observer 
sees that the particle asymptotically approaches R g , but 
never crosses it; the particle distance to the center is 
given by [(5(| 



R(t) = R g + const x exp(—ct/R g 



(34) 



Radius of the collapsing bubble obeys similar asymptotic 
behavior; only the coefficient in the exponent differs by 
a factor of the order of one. 

One should mention that if we neglect gravity the bub- 
ble would collapse into the origin and disintegrate on a 
time scale of a few times of the light crossing time R/c 
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FIG. 7: Radius of the bubble R as a function of time 
t during collapse of initialy static bubble with the radius 
Ro = 941, initial redshift z = 3.07 and a = 100 (solid 
curve). Dashed line shows exponentially decaying fit R(t) = 
865 + 73exp(-£/2523). The unit of length is h/mc, while the 
unit of time is h/mc 2 . 



pl\ . However, if gravity is included this leads to ap- 
pearance of the event horizon which prevents the bubble 
collapse into the origin. The bubble life-time in the pres- 
ence of gravity becomes much larger than R/c. The large 
value of the bubble life time is a consequence of time di- 
lation produced by gravitation field. In the proper co- 
ordinate frame the bubble decays on a short time scale 
~ R/c. However, in the reference frame of a distant ob- 
server (which is relevant to quasar problem) the life time 
is much longer. 

For the thin- wall bubble the total mass-radius relation 
is given by |&fj, l62| 



M 



AttuR 2 



y/1 - {dR/drf 



- 8n 2 u 2 R 3 



(35) 



where r is the interior coordinate time (dr — Ndt). 
When velocity of the bubble surface dR/dr = the ex- 
pression H35fl reduces to our mass-radius relation 12811 for 
n = 1. 



A. Quantum stabilization 

A classical bubble shrinks towards its gravitation ra- 
dius R g and at t 3> R g /c behaves as a black hole. Such 
black hole-like objects can not explain the nature of in- 
trinsically faint quasars. Here we show however that 
quantum effects suppress the collapse and result in ap- 
pearance of long-lived bubbles (not black holes), stable 
on the scale of the Universe age. So far we considered 
the scalar field as classical and Eq. I|35|l defines the clas- 
sical equation of motion of the bubble radius. However, 
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quantum corrections can be crucial for the bubble stabil- 
ity since the later is determined by small terms in Eqs. 
(I10(l - I|12(l . It is known that soliton bubbles which are un- 
stable in the classical model can be stabilized by quantum 
corrections |63l l6t| . From the first glance it is somewhat 
counter-intuitive that for objects 10 Rq in size quantum 
description is appropriate. However, interaction of dark 
matter with incoherent environment is extremely weak. 
Therefore, a coherent state of dark matter can live for a 
long time, which indicates that the quantum description 
is probably adequate. 

To include quantum effects it has been suggested that 
the expression Il35[l for the mass-energy be interpreted as 
the canonical hamiltonian of the bubble at the quantum 
level [H2)EHEB]- bubble wave function ^(R, r) sat- 
isfies the following quan tum mechanical equation in 1+1 
dimensions (h = 1) 62] : 



d 



87rVi? 3 



d 2 

dR 2 



- 167rVi? 4 



V{R,t) = 0. 
(36) 

This equation is interpreted as a one-dimensional wave 
equation for a relativistic particle moving along the semi- 
axis R > with the boundary condition *&(R = 0) = 0. 
Separating variables in Eq. 1361) le ads to the stationary 
one-dimensional wave equation [67j 



{E + 8n 2 u 2 R 3 ) 2 



dR 2 



167rVi? 4 



= 0. (37) 



This equation possesses stationary solutions that are not 
possible in the classical model and the energy spectrum of 
the bubble turns out to be discrete [6(||63- Such station- 
ary quantum bubbles can explain the nature of intrinsi- 
cally faint quasars. At the quantum level the collapse is 
prevented by the uncertainty principle which forbids the 
exact localization of an extended object in space. The 
uncertainty principle yields an effective quantum pres- 
sure which balances the surface tension and gravitational 
attraction producing stationary configurations. The dis- 
crete energy spectrum can be obtained from the quasi- 
classical quantization condition 



Ro 



p(R)dR + ir/4 = (fc + l)7r, 



where momentum p(R) is given by 



(38) 



p(R) = y(E + 8w 2 u 2 R 3 y - 16tt 2 u 2 R\ (39) 

k = 0, 1, 2, ... is the level number and Ro is the classical 
turning point at which p(R) = 0. The turning point is 
determined by the classical mass-radius relation l|28|l . In 
the limit of weak gravity, uR <C 1, Eq. 138|) yields the 
following expression for the energy levels 66] 



£ = 3.59u 1/3 (fc + 3/4) 2/3 . 



(40) 



Bubbles of non-negligible gravitational redshift corre- 
spond to highly excited stationary states with k ~ 



(a/10) 4 ^> 1. For such states the spacing between con- 
secutive levels is small compared to the energy and the 
level distribution can be treated as quasi-continuous. 

Let us now estimate the decay rate of the quantum 
bubble. Partially the picture is analogous to the hydro- 
gen atom where electron moves in the field of nucleus. 
If initially the electron is localized in a small region at a 
large distance R from the nucleus then its further motion 
would be similar to the motion of a classical point par- 
ticle. The originally resting electron will start to move 
radially and fall down to the nucleus during the time 
t c ~ yjR/a, where a is the initial acceleration. How- 
ever, if initially the electron is in a highly excited but 
stationary state its wave function is delocalized. The un- 
certainty principle pressure supports such electron cloud 
from collapse. The decay of the state occurs by means 
of consecutive transitions between stationary states with 
smaller principal quantum numbers and energy loss by 
photon emission. Such decay time t has nothing to do 
with the time t c of the classical radial fall of the electron 
on the nucleus. Nevertheless, t can be estimated from 
classical equations as the time of energy loss by a point 
electron performing radial oscillations around the nucleus 
and emitting electromagnetic waves (even if in the quan- 
tum picture there are no such oscillations). Such a rule 
is a manifestation of the Bohr correspondence principle 
This yields t ~ 4/(R/c) 3 ~ 10" 8 s which agrees 
with the quantum mechanical answer for the decay rate 
of the hydrogen atom excited states, while t c ~ 10~ 16 s. 

The decay of an excited stationary state of the quan- 
tum bubble occurs by means of scalar particle emission. 
We estimate the decay rate using the classical picture as 
the time of energy loss by the classical bubble with the 
radius R(t). We are interested in the bubble life-time 
measured by an outside distant observer. Any correct 
quantum description cannot avoid the Bohr correspon- 
dence principle. So, if the classical bubble treatment, 
mentioned above, yields very long life-time (^> R/c), the 
quantum bubble must also be very long-lived. 

The outside observer cannot witness crossing of the 
horizon R = R g by the bubble surface. As a result, in 
the reference frame of the distant observer the bubble 
wave function ty(R, t) is equal to zero in the inside re- 
gion R < R g [73. Hence, only the classical trajectory 
R(t) between the turning point Rq and the gravitation 
radius R g can contribute to the bubble decay rate. En- 
ergy loss by the bubble surface becomes substantial only 

when R(t) < R^ 3 (see Appendix A below). In our case 
2/3 

R g ^> R and, therefore, the region of intensive energy 
dissipation is not accessible by the classical trajectory. As 
a result, the energy emission is negligible yielding long- 
lived stationary states. In Appendix A we estimate the 
life-time of a quantum bubble as a time of energy loss 
by the classical bubble with the radius R(t) oscillating 
between Rq and R g ; for an order of magnitude estimate 
we take the oscillation period to be Rq / c. The answer is 
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given by Eq. (|A11|1 which in dimension units reads 



zg(z + 2) 6 R (Rq 
(z + l) 12 c ^ I 



(41) 



where I = h/mc is the surface width and z is the bubble 
redshift at R = Rq. For an axion bubble with zq = 0.3, 
Ro = lO 2 i? and I < 0.7 cm Eq. gTJ yields t > 10 18 yrs 
which is much larger then the age of the Universe. 

Apart from stabilization, the quantum effects lead to 
broadening of emission lines. Let us consider an emis- 
sion line with the wavelength Ao produced by atoms in- 
side the bubble. An observer outside the bubble de- 
tects the line at the wavelength A = A (l + z), where 
z = l/yjl - 2M/R - 1 is the redshift of the bubble in- 
terior. At a given bubble mass M the redshift depends 
on the radius R which for the quantum bubble becomes 
uncertain. Its distribution is determined by the square 
of the bubble wave function ^(i?)! 2 . Hence, the emis- 
sion line is detected at the redshift z(R) with a probabil- 
ity proportional to ^(i?)! 2 and, instead of a sharp line, 
the observer would detect a broaden peak centered near 
z(R) corresponding to a maximum of ^(i?)! 2 . In the 
quasiclassical (WKB) approximation the wave function 
in the interior coordinate frame can be written as follows 

mini 



c 



: sin 




p{x)dx + - 



(42) 



is a solution to the left of the classical turning point Rq. 
In the limit of weak gravity (uR < 1) the normalization 
constant is C ~ 1.24^/ E j Rq. In the classical forbidden 
region R > Rq the wave function decreases as 



9(R) 



(43) 



Near the turning point p(R) oc Rq — R. For R > Rq 7 
\l/(i?) decays fast on the scale AR ~ 1/a 2 / 3 , hence, in 
the classical forbidden region ^f(R) ~ 0. For R < Ro, 
^(R) undergoes fast oscillations. Averaging over the os- 
cillations yields in the vicinity of the turning point 



\m)\ 



2 Q(Ro- R) 



\/ Rq — R 



(44) 



where Q(x) is the step function. The wave function 1441) 
has a peak at R = Rq and, therefore, the external ob- 
server would detect a peak in the bubble radiation spec- 
trum at the wavelength A' = Ao(l + z(Rq)) as if the 
emission line is redshifted according to the classical mass- 
radius relation (J2SJ. Hence, for the quantum bubbles the 
same classical Eqs. (|28|l . (|30|l determine the redshift of 
the emission features. However, the detected emission 
line profile F(X) would be different from those emitted 
by the atom. If for the atom F(X) cx <5(A — Ao) then the 
external observer would detect 



F(X) 



6(A - A') 
a/A - A 7 ' 



(45) 



The line profile becomes asymmetric. However, the line 
width due to quantum broadening is negligible and other 
broadening mechanisms, e.g., atom motion, would prob- 
ably wash out the quantum profile. 

Finally we emphasize the dramatic difference between 
the classical and quantum descriptions. In the classical 
picture the bubble surface localizes after some time at 
the event horizon which leads to black hole formation. 
However, due to the uncertainty principle, quantum me- 
chanics does not allow localization of the bubble surface. 
Eq. H34|) suggests that in the reference frame of a distant 
observer the bubble wave function in the WKB approxi- 
mation at R 



R g behaves as 

|*OR)| 2 oc 0(i * -jRs) 



R — R n 



(46) 



The WKB approximation becomes invalid in a small 
vicinity of R g and the exact wave function remains finite. 
The normalization integral J \^(R)\ 2 dR of the quasiclas- 
sical wave function logarithmically diverges at R — ► R g . 
However, the very slow (logarithmic) divergence suggests 
that if we cut off the integration at any reasonable dis- 
tance to R g the total contribution from the vicinity of 
R g is not very large. Hence, in the quantum description 
there is quite large probability that the bubble radius 
is substantially greater than R g . So, in contrast to the 
classical picture, the long-lived quantum bubble does not 
behave as a black hole. Light escapes from the bubble 
interior, but it possesses gravitation redshift determined 
by Eqs. 1(28]). lj3UJ|. This is our model of axionic quasars. 

One might expect that the 1/(R—R g ) divergence of the 
WKB wave function (|46|l should produce features in the 
bubble radiation spectrum at low frequencies. However, 
the redshift z reduces the radiation power by a factor 
of 1/(1 + z) 2 = (1 - Rg/R) 72]. This compensates the 
wave function divergence and yields no emission features 
at low frequencies. 



IV. THE NATURE OF BRIGHT QUASARS 

The results obtained show that some fraction of 
quasars is probably bubbles of scalar field with periodic 
interaction potential. They are born in active galaxies 
and ejected into surrounding space. The ejection sug- 
gests that the bubble mass must be much smaller then 
the mass of the parent galaxy. A bubble with a mass 
10 9 M Q and gravitational redshift z = 1 has a radius 
R = 10~ 4 pc and, if located at a distance 1 Mpc, would 
have an angular size of 0.00004" which is well below the 
resolution limit of current telescopes in the optical band 
(0.1"). Therefore, the dark matter bubbles can explain 
only the faint quasars with unresolved structure. Al- 
lowed values of the mass and decay constant / of axion 
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fit well into this picture suggesting the axion as primary 
candidate for composition of the dark matter bubbles. 

The question whether or not all dark matter in the 
Universe consists of axions requires further detailed anal- 
ysis. However, existence of bright quasars with resolved 
"host" galaxies suggests that some part of dark matter is 
composed of other particles. Fuzz around bright quasars 
3C 48 (z = 0.37) [zl, 3C 273 (z = 0.158) 3C 279 
(z = 0.536) [75| was found with the redshift equal to the 
quasar redshift. The angular size of the fuzz is 0.5" (3C 
279), 4" (3C 48) and 10" (3C 273). Due to large angular 
size and brightness such quasars can not be axionic bub- 
bles. However, based on observations, Arp has suggested 
that 3C 273 and 3C 279 are members of the Virgo cluster 
21 Mpc away), while 3C 48 belongs to the Local Group 
6j. There is considerable evidence from the morphology 
7g and from X-ray emission [TtJ that this may be the 
case and, hence, those objects possess noncosmological 
redshift. 

Dark matter bubbles or droplets with radii 10 — 500 
pc can explain the observed angular size of the bright 
quasars provided that the redshift is gravitational. How- 
ever, if we estimate the required mass from Eq. I|30|) 
we obtain 10 14 — 10 15 M Q which is much greater than a 
typical mass of a large galaxy 10 12 M©. Such big masses 
definitely would not fit into any realistic picture of galaxy 
clusters and contradict to the scenario of quasar ejection 
from galaxies. 

Here we discuss a solution of the Einstein-Klein- 
Gordon equations describing an object with galactic 
size, large gravitational redshift, but possessing negligi- 
bly small mass compared to the mass of a galaxy. The 
solution is similar to a spherical capacitor known in elec- 
trostatic. If one plate of the "capacitor" has the mass 
+M, while the other one possesses the mass —M the total 
mass and gravitational field outside the capacitor is equal 
to zero. However, the capacitor interior has a nonzero 
gravitational potential and, hence, nonzero gravitational 
redshift. Such a capacitor can be realized for scalar fields 
with negative (attractive) interaction potential. The field 
kinetic energy provides positive contribution to the en- 
ergy density, while the contribution from interaction is 
negative. Space regions where the kinetic energy domi- 
nates play a role of the plates with positive mass, while 
interaction dominating regions are analogous to the neg- 
ative plates. One can anticipate the effect from the struc- 
ture of the Einstein equation i|ll|) for the metric compo- 
nent N that determines the space redshift z = 1/N — 1. 
If V < then the combination ip' 2 — 2g 2 V everywhere 
gives the positive contribution to N producing large to- 
tal redshift. However, according to Eq. (|16f) . the mass 
density depends on ip' 2 + 2g 2 V and if V < the total 
mass can be very small. 

To demonstrate the effect we consider a complex scalar 
field ip with the simplest attractive potential 



V(ip) = -m 2 \ip\ 2 /2. 



(47) 



with the speed greater then the speed of light c, possess 
momentum p > mc and the energy E = \Jp 2 c 2 — m 2 c 4 . 
There is an evidence that neutrino could be a tachyon 
[78l |7jJ , which makes the hypothesis that dark matter is 
partially composed of tachyons not so exotic. Tachyons 
as a candidate for dark matter and dark energy are cur- 
rently discussed in the literature (see, e.g., 0, 
and references therein). They can explain the observa- 
tions at both large and galactic scales [8^, Q • 

The equilibrium field configurations are those in which 
the metric is time independent. The scalar field ip itself 
can oscillate with frequency u>, ip(t, r) = ip(r) exp(— iut), 
however due to U(l) symmetry of the Lagrangian the 
space-time geometry is static. In dimensionless units, 
defined by Eq. iJSJ , the stationary Klein-Gordon and Ein- 
stein equations for ip(r) are 0] 



if)' (g 2 + 1 

7^2 



Ip" UJ 2 1p 



m 2 N 2 



+ V = 0, (48) 



N' = 



9 



N 



9 2 -l 



+ r [ 4>' 2 - g' z i/j z + 



2„;,2 



Lo 2 g 2 ip 2 



2 N 2 



, (49) 



(50) 



For simplicity we assume ui — 0. In the limit of zero 
gravity the Klein-Gordon equation 



Alp + ip = 
has a spherically symmetric solution 



i> = ipo 



smr 



(51) 



(52) 



that describes a droplet with specially oscillating energy 
density E(r) which, at r ^> 1, is given by E{r) — ip /2 /2 — 
i/j 2 /2 — -0q cos(2r)/2r 2 . The solution is analogous to a 
set of concentric spherical plates charged alternatively 
with positive and negative mass. 

General relativistic treatment modifies the formula. 
First we discuss the solution of Eqs. I|48fl - i|50|) in the 
limit \tp\ <C 1. In such limit one can take g sa 1, then Eq. 
gSJ leads to 



ip' ( - + rip 2 + ip" +ip = Q. 



(53) 



Further we take the nonlinear term rip 2 in the zero order 
approximation (|52|l . This term is important at r 3> 1 and 
yields a decay of the scalar field faster then 1/r. Also we 
average the term rip 2 over the period of space oscillation, 
that is assume rip 2 w ipQ/2r, and finally obtain 



(2 + ipl/2)iP'/r + ip" + i/j = 0. 



(54) 



Solution of this equation is expressed in terms of the 
Bcsscl function 



The free field describes tachyons that always propagate 



iP 



^^0^1/2+02/4 (r) 



(55) 
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FIG. 8: Redshift z(r) of space as a function of distance r to 
the droplet center. Note that the redshift is normalized to 
its center value z c and the scale along the horizontal axis is 
logarithmic. The unit of length is h/mc. 



and has the following asymptotics: ip m ipo at r <C 1, and 
ip x, ip smr/r 1+ ^o/ 4 at r 3> 1. The total mass-energy of 
the droplet is (in units m^/m) 



E = 



r /)2 

9 



dr < O 4 , 



(56) 



which for ip <C 1 is negligible. However, for the gravita- 
tional redshift at the droplet center z c we obtain 



exp 



r(ip' 2 + ip 2 )dr 



1 = e - 1 = 1.718, 



the value is independent of ipo (in the limit ip$ <C 1). 
If we move away from the droplet center the redshift of 
the space decreases from z c to zero with the asymptotic 
behavior z(r) = l/r^'o/ 2 (r > 1). 

Numerical solution of Eqs. Q48 [) -Q5U [I shows that de- 
pending on the scalar field at the droplet center ipo the 
gravitational redshift can have any large values. In Fig. 
|S]we plot the spatial distribution of the redshift z(r) for 
different values of ipo. Near the droplet center the red- 
shift changes on the scale of a few H/mc, while far from 
the center the spatial variation becomes very weak. 

Our model can provide an estimate of the character- 
istic scales. If, e.g., we take m — 0.64 x 10~ 23 eV 
then the unit of length is h/mc = 1 pc. A bright 
baryonic nucleus located a distance r = 10 kpc from 
the center of the tachyonic droplet with ipo = 0.2 (the 
droplet energy \E\ < ip^m^/m = 10 10 M Q ) would pos- 
sess a redshift z — l/r^°/ 2 = 0.83. Fuzz surround- 
ing the baryonic nucleus with a diameter Ar = 200 
pc would have a dispersion of the gravitational redshift 
Az w VoAr/2r 1+ ^o/ 2 = 0.00033 = 100 km/s, the red- 
shift increases towards the droplet center. Such system- 
atic change in the redshift across the fuzz can be directly 



measured and serve as a test of our theory. However, the 
position of the bright nucleus might also coincide with 
the droplet center, such a picture can be realized pro- 
vided we choose for m a smaller value. 

Our finding shows that substantial gravitational red- 
shift of a large space volume does not necessarily require 
presence of a big mass. Even a volume of galactic scale 
can possess any value of the gravitational redshift pro- 
duced by dark matter with the total mass negligibly small 
compared to galactic mass. Ordinary baryonic matter 
placed in such a volume behaves as being possessed of 
an intrinsic redshift. The effect can explain the nature 
of intrinsically bright quasars as nuclei of forming nearby 
small galaxies embedded in a clot of scalar field with neg- 
ative interaction. It also can account for recent observa- 
tion of the field surrounding the Seyfert galaxy NGC 7603 
where four galaxies with substantially different redshifts 
are apparently connected by a narrow filament |85| . 

We want to emphasize that the interaction potential 
(P7|) provides a simple demonstration of the effect, rather 
then a precise quantitative description of the phenomena. 
Realistic potential must be derived based on detailed 
study of bright quasars including their environment and, 
probably, should give rise to a faster asymptotic decay of 
the scalar field with the distance from the droplet center. 



V. DISCUSSION 
A. Axions 

In this paper we argue that the problems of quasars 
and dark matter are mutually related. We found that 
bubbles of scalar field with periodic interaction poten- 
tial can explain the nature of intrinsically faint point-like 
quasars associated with nearby galaxies. Typical abso- 
lute magnitude M v of such quasars lies in the range —8 to 
-13 (optical luminosity L — 10 5 — 1O 7 L ) (Hill, the 
redshift is mostly gravitational and found to be quan- 
tized. The bubbles are born in nuclei of active galaxies 
and ejected into surrounding space. They cluster at a dis- 
tance upto 100—300 kpc from the parent galaxy. About 
15 such objects have been discovered close to M82, the 
nearest active galaxy to the Milky Way 01 1 an d about 
10 in the vicinity of NGC 3628 ||. 

The observed five peaks in the quasar redshift distri- 
bution match well the theoretical result with only one 
free parameter, which is a strong argument in favor of 
our theory. The hypothetical axions, one of the leading 
dark matter candidate, fit well into the quasar picture 
and can account for the bubble composition. The axion 
mass range constrained by astrophysical and cosmologi- 
cal arguments yields the necessary bubble life-time and 
the size which agrees with the observed quasar bright- 
ness. 

Properties of the intrinsically faint QSOs, com- 
bined with equations for the bubble mass M = 
0.00752a 2 m 2 1 /m = 2.94m(eV) x lO n M , radius R = 
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FIG. 9: Astrophysical and cosmological exclusion regions 
(hatched) for the axion mass m a . The dotted "inclusion re- 
gion" indicate where axions could plausibly be the cosmic 
dark matter. Axion mass determined from quasar observa- 
tions fits in the open window. Also shown is the projected sen- 
sitivity range of the search laboratory experiments for galactic 
dark matter axions. 



0m29a 2 h/mc = 2.73m(eV) x 10 6 i? Q and the relation 
© m = 0.62 eVxl0 7 GeV// [Hi, allow us to deter- 
mine the axion mass m. The quasar luminosity suggests 
that the bubble radius is larger then 1O 3 J?0 which yields 
in > 4 x l(T 4 eV and M > 10 8 M Q . From the other hand, 
the quasar ejection from active galaxies implies that the 
bubble mass M must be much smaller then the galactic 
mass. It is reasonable to constrain M < 10 9 M Q which 
leads to m < 3 x 10~ 3 eV and R < 10 4 i? Q . We conclude, 
the axion mass is m = 0.4 — 3 meV. This value fits in 
the open window for the axion mass constrained by as- 
trophysical and cosmological arguments, as displayed in 
Fig. [5J which unambiguously points towards the axionic 
nature of dark matter composing the intrinsically faint 
point-like quasars. One can see that current cavity ax- 
ion search experiments in Livermore [Sfil l87l ] and Kyoto 
University |30J are looking for the axion in an unlikely 
mass range which deviates by two orders of magnitude 
from our result (see Fig. EJ. Probably now, when the 
axion mass is established from quasar observations, the 
axion has a better chance to be discovered. One should 
mention that radio telescopes are suited to search for ax- 
ions of higher mass and allow searches in the range 
obtained in our paper. 

Understanding the mechanisms of quasar formation re- 
quires further detailed study. Bose-Einstein condensa- 
tion of axions in galactic centers can be a possible natu- 
ral mechanism of bubble creation. Dark matter axions, if 
they exist, form halos around galaxies. The halo axions 
are in a quantum degenerate non-equilibrium regime with 
the de Broglie wavelength ~ 1 m which is much larger 
then the interparticle spacing [23] • The ground state of 
the system would be a collapsed axion condensate at the 
galactic center. However, the kinetics of condensate nu- 
cleation is governed by two factors: inter-particle colli- 
sions and occupation number of the condensate particles. 
Due to extremely weak interactions between axions the 



collision rate is very small, so the galactic halo remains 
metastable for a long time. 

Recent experiments on Bose condensate formation in 
magnetic traps [8fl| show that in a weak cooling regime 
the condensate nucleation occurs in two stages: slow lin- 
ear growth which then, at a critical point, triggers fast 
exponentially growing instability, caused by the effect of 
bosonic stimulation. Similar situation might take place in 
nuclei of active galaxies. At the first stage, the axions are 
slowly (during million years) accumulated at the central 
part of the galaxy. When the mass of such a ball becomes 
critical (particle occupation number is big enough), it 
triggers exponential instability due to bosonic stimula- 
tion. A dense condensate cloud, coherent on astronomi- 
cally large scale, starts to form rapidly. The cloud then 
collapses under its own gravity. Such a mechanism is 
similar to type la supernovae, when a star explodes af- 
ter accumulation a critical mass, and suggests formation 
of objects with approximately equal masses. Why does 
the cloud collapse lead to formation of bubbles? The 
answer comes from a three-dimensional numerical sim- 
ulation of the evolution of inhomogencities in the axion 
field. Such a simulation has indeed demonstrated forma- 
tion of bubble-like structures (see Fig. 5a in Ref. [Ofij). 
Gravitational cooling is probably an important processes 
involved in quasar nucleation |9l| . 

Finally we want to mention a new (preliminary) ex- 
perimental result which became available after our the- 
ory was first presented in DARK 2004 [H|. The PVLAS 
experiment on axion laser production and detection ap- 
parently observes a sig nal which might be produced by 
axion- like particles 92] . If the observed signal is indeed 
caused by axions then the particle mass measured in the 
PVLAS experiment is m = 1.0 ± 0.1 meV 92]. This 
value is right in the middle of the axion mass interval 
0.4 — 3 meV we predicted based on quasar observations! 
One should mention that another parameter, the axion- 
photon coupling, measured by PVLAS is irrelevant to 
our theory. Our estimate of the axion mass is based 
only on the relation J2J) between the global symmetry- 
breaking scale / and m. From this perspective one can 
treat our and PVLAS results as mutually complemen- 
tary. Only combination of the two independent results 
confirm, based on observations, the relation J5J which 
is the key ingredient of any axion model. This provides 
an evidence that the new particle, detected by PVLAS 
and responsible for the quasar redshift quantization, is 
indeed axion and not another yet unknown pseudoscalar 
particle. 



B. Tachyons 

Existence of the both intrinsically faint point-like and 
bright resolved quasars associated with galaxies suggests 
that dark matter in the Universe is probably composed 
of several species. We found that tachyonic clots can 
produce substantial gravitational redshift in a galactic 
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scale even when the total mass of the object is negligi- 
ble. Bose condensate of tachyons can explain the nature 
of intrinsically bright quasars as forming galactic nuclei 
emerged into droplets of scalar field with negative inter- 
action. The droplets possess large gravitational redshift 
on kpc scales. About 42 such objects have been discov- 
ered in the vicinity of the Seyfert 1 galaxy NGC 6212 
which is approximately 120 Mpc away |93| . They form 
a cloud of QSOs with a characteristic size 2 Mpc and 
absolute QSO magnitudes from —20 to —15. 

Being large and luminous, the tachyonic objects domi- 
nate in samples of quasars associated with relatively dis- 
tant, > 50 Mpc, active galaxies. At such distances the ax- 
ionic bubbles are so faint that they are rarely detectable 
as individual sources. According to our theory, the red- 
shift quantization is the property of axionic bubbles and 
must not be present for bright tachyonic quasars. In 
agreement with the theory, no periodicity in the intrin- 
sic redshift distribution has been found in a sample of 
tachyon dominating quasars associated with the distant 
(0.01 < z ga i < 0.3) galaxies H3- 

If within a distance d < 10 Mpc there are, at least, 
several bright tachyonic quasars, then inside the region 
d < 100 Mpc thousands of them must be detected. Prob- 
ably such objects dominate in the recent First Bright 
Quasar Survey which radio-selected sources brighter than 
18 optical magnitude [13 HU; axionic bubbles are usu- 
ally too faint and not selected in such a survey. Luminous 
tachyonic quasars can be also seen at cosmological dis- 
tances where the cosmological redshift contribution be- 
comes substantial. The gravitational redshift of tachy- 
onic quasars can also be very small which provides a con- 
tinuous connection between quasars and active galactic 
nuclei, the classical quasar explanation. In general case, 
however, the quasar redshift possesses substantial gravi- 
tational component and, hence, is not a measure of the 
object distance. 

Proximity of some quasars and galaxies with approxi- 
mately equal redshift indicates only that for some objects 
the gravitational redshift is small compared to cosmolog- 
ical. Meanwhile, the statistical evidence, mentioned in 
the Introduction, implies that the gravitational redshift 
dominates for the majority of quasars in the selected sam- 
ple. It is also worth to note that if a dwarf galaxy lies at 
a distance of a few kpc from the tachyonic quasar, both 
such objects are emerged into the same tachyonic clot 
and, as a consequence, must have close redshifts. This 
could be a reason why some quasars and associated galax- 
ies possess equal redshifts and account for the nature of 
binary quasars. 

It is known that quasars violate the Hubble redshift- 
apparent magnitude relation which is an argument 
against the cosmological redshift. The argument would 
be very strong if the quasar apparent magnitude m v 
proved to be a good indicator of their distance. If objects 
with equal luminosity possess a uniform space distribu- 
tion then their number N detected per unit m v obeys the 
law \og(dN/dm v ) — const + 0.6m v j9|J. In Fig. ITU1 we 
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FIG. 10: Logarithm of the number of quasars within the 
Am, = 0.5 magnitude bins in the First Bright Quasar Sur- 
vey of the north Galactic cap as a function of the apparent 
magnitude m v (dots) ;94]. The solid line is the best linear fit. 



plot \og(dN / dm v ) as a function of the optical magnitude 
Tfi v for quasars discovered by the First Bright Quasar 
Survey (dots) 94]. The solid line is the best linear fit 
which yields the slope of 0.58 very close to 0.6. Such an 
agreement suggests that for a large sample mostly the 
distance determines m v of the bright quasars and, hence, 
the scattered Hubble z — m v diagram implies noncosmo- 
logical redshift for the majority of quasars in the Survey. 

Arp has noticed an empirical sequence of quasar devel- 
opment in which initially point-like objects at relatively 
faint apparent magnitudes transform into lower-redshift 
compact objects with "fuzz" around their perimeters, 
and then into small, high surface brightness galaxies 
P. loH losT l . In the light of our theory this means that faint 
axionic bubbles, which are born in and ejected from ac- 
tive galaxies, later serve as nucleation centers for tachyon 
Bose-Einstein condensation. Condensate formation leads 
to appearance of a large in size, bright quasar which 
then gradually transforms into a small companion galaxy. 
Tachyonic clots can account for the systematic excess of 
the redshift of small companion galaxies compared to the 
redshift of the dominant galaxy [|| . Continuity suggests 
that the dark halos of galaxies, known from their rota- 
tion curves, must be remnants of the tachyonic matter. 
Possibility of such scenario is a subject of future study. 

The process of axionic quasar nucleation could produce 
substantial disturbances in galaxies. Such catastrophic 
events can contribute to the observed gamma-ray bursts. 
However, existence of nearby quasars suggests that such 
events also occur in our epoch, not only in the early Uni- 
verse. If we divide the total number of galaxies in the 
Universe, 8 x 10 10 , by the frequency of the gamma-ray 
bursts (~ 1 per day) we obtain that in average the catas- 
trophic event occurs with the interval of 100 million years 
per galaxy. This value is of the order of the time span be- 
tween geological periods on Earth which might indicate 
on their Galactic origin. More often catastrophic events 
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occurring in the Local galactic group might be responsi- Mars, can verify this hypothesis. 

ble for the change in the geological ages within a period. I wish to thank A. Belyanin and A. Sokolov for helpful 

Study of the geological periods on other planets, e.g., on discussions. 
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APPENDIX A: ENERGY EMISSION FROM A SHRINKING BUBBLE 

Here we calculate energy loss by a shrinking spherically symmetric bubble caused by emission of scalar particles. 
For an order of magnitude estimate one can omit the effect of gravity. Then the evolution of the scalar field ip(t, r) is 
described by sine-Gordon equation 



Cp — ip + sin ip = 2(p I r, 

where r is the radial coordinate. Without right-hand side, Eq. I|A1|) has an exact, so-called kink, solution 



ipo = 4 arctan < exp 



± 



(r-vt- Ro) 



VY^v 2 



(Al) 



(A2) 



where Rq 3> 1 is the initial bubble radius. The solution describes a kink (space region where if changes from 2ir to 0) 
propagating with constant velocity v; the kink's size is I ~ \/l — v 2 . 

If I <C R(t), where R(t) is the bubble radius, r.h.s. of (|A1|) may be treated as a small perturbation. Eq. (|A1|) 
possesses approximate solution in the form of the kink l|A2f) with parameters slowly changing in time under the action 
of the perturbation. In particular, the kink shrinks due to its surface tension so that the bubble radius and the 
velocity evolve as J9!§ 



R(t) = R cn(V2t/R a , 1/V2), v(t) = Jl- R A {t)/Rt 



(A3) 



where cn stands for the elliptic cosine with the modulus 1 / \[2. Such a process is accompanied by emission of scala r 
particles which yields the energy loss. We estimate the energy loss following the original work of Malomed [o^. llOCj . 
In terms of the inverse scattering technique, the spectral density of the emitted energy E e (t, q) is 



(A4) 



where q is the radiati on w avenumber and the perturbation-induced evolution equation for the complex amplitude 
B(t.a) is given bv [99llTooj 



dB 
~dt 



2(A 2 + 7 2 ) J_ 



dr [X 2 -j 2 ~ 2iA7tanh 



r — vt 



exp + q 2 t — iqr^J d r ip 0l (A5) 



where A = y/l + q 2 — q and 7 = (1 + v)/2y/l — v 2 . Calculating the integral in 1A5(I yields 

dB _ in [A 2 (l - «)(1 - vTTu) - v/2] ex P (V 1 + 9 2 * " l< l vt ) 
~dt ~ (l + t;)/4 + A 2 (l-w) cosh [irqVT^ 2 '/2] 

If v slowly varies with time one can take v » const in Eq. I|A6|I . then after integration we obtain 

7T [a 2 (i -v)(i- VTTv) - v/2] ex P (W 1 + <i 2t - iqvt) - 1 



(A6) 



Therefore 



[(1 + u)/4+ A 2 (l - v)] (y/1 + q 2 - qv) cosh [vrgVl - v 2 /2] 
dE e _ 167r[A 2 (l-w)(l-VTT^)-«/2] 



(A7) 



' 2 sin 2 



y/1 + q 2 - qv) t/2 



dq [(l + v)/4 + X 2 {l-v)] 2 (^/TTq 2 -qv) 2 cosh 2 [ 7 rgVT 3 ^ 2 /2] 



(A8) 
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Integration of 1|A8(1 over dq gives the emitted energy as a function of time E e (t) — dq(dE e J dq) . In Eq. (|A8fl sine 

is a fast oscillating function, so we substitute sin 2 (a;) — ► 1/2. The radiation power increases when the kink's velocity 
v approaches the speed of light c = 1. Assuming 1 — v <C 1, integration of Eq. (|A8|) yields 

« ! >*(d!W "™(m) ° (A9) 

The emitted energy becomes comparable with the initial bubble energy E n = 8Rq when the bubble radius reaches 
the value i?* « R^ 3 - This value agrees with those obtained in [6l|. 

If the bubble shrinks from Rq to the gravitation radius R g S> R* the radiated energy is 



R \\ (zq + 1) 12 £q 
J z 6 (z + 2) 6 i? 2: 



7.10(^) « ,^„ -£, (A10) 



where z = 1/ yT~ Rg/Ro — 1 is the initial bubble redshift. To emitt all its energy the bubble must oscillate between 
R g and Rq about Eq/E 9 cycles. As a result, the bubble life-time is 



E Q zg(z + 2) 6 flg 

^ ~ (zo + iy* ' (A11) 



which for z = 0.3 yields i - 0.005i?g. 
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